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Abstract 



We present the first direct measurements of the pion valence quark momen- 
tum distribution which is related to the square of the pion light-cone wave 
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function. The measurements were carried out using data on diffractive dis- 
sociation of 500 GeV/c 7r~ into di-jets from a platinum target at Fermilab 
experiment E791. The results show that the \qq) light-cone asymptotic wave 
function, which was developed using perturbative QCD methods, describes 
the data well for Q 2 ~ 10 (GeV/c) 2 or more. We also measured the trans- 
verse momentum distribution of the diffractive di-jets. 
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The internal momentum distributions of valence quarks in hadrons enter the calcula- 
tion of a large variety of processes such as electroweak decays, diffractive processes, meson 
production in e + e _ and 77 annihilation, relativistic heavy ion collisions, and many others 
[1]. The momentum distribution amplitudes are generated from the valence light-cone wave 
functions integrated over k t < Q 2 , where k t is the intrinsic transverse momentum of the 
valence constituents and Q 2 is the total momentum transfer squared (Eqn. 5). Because of 
the close relationship between the two, the distribution amplitudes are often referred to as 
the light-cone wave functions [2]. Even though these amplitudes were calculated about 20 
years ago, there have been no direct measurements until those reported here. Observables 
which are related to these distributions, such as the pion electromagnetic form factors, are 
rather insensitive to the light-cone wave functions. 

The pion wave function can be expanded in terms of Fock states: 

* = a^qq) + a 2 \qqg) + a 3 \qqgg) + ■■■ . (1) 

The first (valence) component is dominant at large Q 2 . The other terms are suppressed by 
powers of l/Q 2 for each additional parton, according to counting rules [2,3]. In contrast, 
parton distribution functions are inclusive momentum distributions of partons in all Fock 
states. Here we are concerned with the momentum distribution of only the valence quark- 
antiquark part. 

Two functions have been proposed to describe the momentum distribution amplitude for 
the quark and antiquark in the \qq) configuration. The asymptotic function was calculated 
using perturbative QCD (pQCD) methods [4-6], and is the solution to the pQCD evolution 
equation for very large Q 2 (Q 2 — > 00): 

<j> aa (x)=V^x(l-x). (2) 

x is the fraction of the longitudinal momentum of the pion carried by the quark in the 
infinite momentum frame. The antiquark carries a fraction (1 — x). Using QCD sum rules, 
Chernyak and Zhitnitsky (CZ) proposed [7] a function that is expected to be correct for low 
Q 2 : 

<p cz (x) = 5V3x(l-x)(l-2x) 2 . (3) 

As can be seen from Eqns. 2 and 3, and from Fig. 2, there is a large difference be- 
tween the two functions. Measurements of the electromagnetic form factors of the pion were 
considered to be the best way to study these wave functions. A comprehensive summary 
of the status of these measurements was published recently [2]. Both existing methods of 
measuring the pion electric form factor suffer from major drawbacks. Those done using 
pions elastically scattered from atomic electrons measure the pion electric form factor only 
for very low Q 2 . The alternative method is to measure the electron-pion quasi- free scat- 
tering cross section [2]. The relation between the form factor and the longitudinal cross 
section is model-dependent as it includes the p — > nii + matrix element. Finally, the form 
factor is related to the integral over the wave function and the scattering matrix element, 
reducing the sensitivity to the wave function. Indeed, as shown in [2] both wave functions 
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can be made to agree with the experimental data. A similar situation exists for inelastic 
form factors, decay modes of heavy mesons, etc. Recent model-dependent analyses of CLEO 
data on meson-photon transition form factors [8,9] are consistent with the asymptotic wave 
function. The problem is that comparisons to these observables are not sensitive to details 
of the wave function and thus cannot provide critical tests of their x- dependence. Another 
open question is what can be considered to be high enough Q 2 to qualify for perturbative 
QCD calculations, what is low enough to qualify for a treatment based on QCD sum rules 
and how to handle the evolution from low to high Q 2 . 

In this work we describe an experimental study that maps the momentum distribution 
of the q and q in the \qq) Fock state of the pion. This provides the first direct measurement 
of the pion light-cone wave function (squared). The concept of the measurement is the 
following: a high energy pion dissociates diffractively on a heavy nuclear target imparting 
no energy to the target so that it does not break up. This is a coherent process in which 
the quark and antiquark in the pion break apart and hadronize into two jets. If in this 
fragmentation process each quark's momentum is transferred to a jet, measurement of the 
jet momenta gives the quark and antiquark momenta. Thus: 



The diffractive dissociation of high momentum pions into two jets can be described, like 
the inclusive Deep Inelastic Scattering (DIS) and exclusive vector meson production in DIS, 
by factoring out the perturbative high momentum transfer process from the soft nonper- 
turbative part [10]. This factorization allows the use of common parameters to describe 
the three processes. The virtuality of the process, Q 2 , is given by the mass-squared of the 
virtual photon in inclusive deep inelastic scattering (DIS) and in exclusive light vector me- 
son production in DIS. In exclusive DIS production of heavy vector mesons this mass is 
proportional to that of the produced meson (e.g., J/ijS). For diffractive dissociation into two 
jets, the mass-squared of the di-jets plays this role. From simple kinematics and assuming 
that the masses of the jets are small compared with the mass of the di-jets, the virtuality 
and mass-squared of the di-jets are given by: 



where kt is the transverse momentum of each jet and reflects the intrinsic transverse momen- 
tum of the valence quark or antiquark. By studying the momentum distribution for various 
kt bins, one can observe changes in the apparent fractions of asymptotic and Chernyak- 
Zhitnitsky contributions to the pion wave function. 

Fermilab experiment E791 [11] recorded 2 x 10 10 events from interactions of a 500 GeV/c 
7r~ beam with carbon (C) and platinum (Pt) targets. The trigger included a loose require- 
ment on transverse energy deposited in the calorimeters. Precision vertex and tracking 
information was provided by 23 silicon microstrip detectors (6 upstream and 17 downstream 
of the targets), ten proportional wire-chamber planes (8 upstream and 2 downstream of the 
targets), and 35 drift-chamber planes. Momentum was measured using two dipole magnets. 
Two multicell, threshold Cerenkov counters were used for ir, K, and p identification (not 



•^measured 



Pjetl 



(4) 



Pjetl + Pjet2 



Q 2 ~M 2 = 




(5) 



A 



needed for this analysis). Only about 10% of the E791 data was used for the analysis pre- 
sented here. From these data, we selected interactions which were exclusive two jet events. 
This focused on the jets as materializations of the valence quark and antiquark in the pion. 

The data were analysed by selecting events in which 90% of the beam momentum was 
carried by charged particles. This reduced the effects of the unobserved neutral particles 
and allowed for precise measurement of transverse momentum. The selected events were 
subjected to the JADE jet-finding algorithm [12]. The algorithm uses a cut-off parameter 
( m cut) whose value was optimized for this analysis using Monte Carlo simulation studies 
in order to optimize the identification of di-jets. The di-jet invariant mass was calculated 
assuming that all the particles were pions. To insure clean selection of two-jet events, a 
minimum k t of 1.25 GeV/c was required. Furthermore, the di-jet nature of these events 
was verified by examining their relative azimuthal angle, which for pure di-jets should be 
180°. Strong peaking at 180° was observed (FWHM ~ 5°), and only events within 20° of 
back-to-back were accepted for this analysis. 
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FIGURES 
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FIG. 1. q% distributions of di-jets with 1.25 GeV/c < k t from interaction of 500 GeV/c tt~ 
with Carbon and Platinum targets. 

Diffractive di-jets were identified through the e~ bg * dependence of their yield (qf is the 
square of the transverse momentum transferred to the nucleus and b = <R 3 > where R is the 
nuclear radius). Fig. 1 shows the q\ distributions of di-jet events from platinum and carbon. 
The different slopes in the low q\ coherent region reflect the different nuclear radii. Events 
in this region come from diffractive dissociation of the pion. 

The basic assumption that the momentum carried by the dissociating qq is transferred to 
the di-jets was examined by Monte Carlo (MC) simulation. MC samples with 4 and 6 GeV/c 2 
mass di-jets were generated with two different x dependences at the quark level. The x- 
distributions were calculated by squaring the asymptotic and the Chernyak-Zhitnitski (CZ) 
wave functions. One sample was simulated using the asymptotic wave function and the other, 
the CZ function. The four samples were allowed to hadronize using the LUND PYTHIA- 
JETSET [13] package and then passed through a simulation of the experimental apparatus 
to account for the effect of unmeasured neutrals and other experimental distortions. 
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FIG. 2. Monte Carlo simulations of squares of the two wave functions at the quark level (left) 
and of the reconstructed distributions of di-jets as detected (right). 4>\ sy is the asymptotic function 
(squared) and 4>q Z is the Chernyak-Zhitnitsky function (squared). The di-jet mass used in the 
simulation is 6 GeV/c 2 and the plots are for 1.5 GeV/c < k t < 2.5 GeV/c. 

In Fig. 2 the initial distributions at the quark level are compared with the final dis- 
tributions of the detected di-jets, including distortions in the hadronization process and 
influence due to experimental acceptance. As can be seen, the qualitative features of the 
two distributions are retained. The results of this analysis come from comparing the ob- 
served x-distribution to a combination of the distributions shown, as examples, on the right 
of Fig. 2. 

For all results in this paper, we used data from the platinum target as it has a sharp 
diffractive distribution and a relatively low background. It is also expected that due to the 
color transparency effect [6,14,15] this heavy target will better filter out the high Fock states. 
We used events with g 2 < 0.015 GeV/c 2 . For these events, the value of x was computed from 
the measured longitudinal momentum of each jet (Eqn. 4). A background, estimated from 
the x distribution for events with larger qf, was subtracted. This analysis was carried out in 
two windows of k t : 1.25 GeV/c < k t < 1.5 GeV/c and 1.5 GeV/c < k t < 2.5 GeV/c. 
The experiment data were compared to Monte Carlo simulations of di-jets having a mass 
of 4 GeV/c 2 for the lower window and 6 GeV/c 2 for the higher window. The resulting x 
distributions are shown in Fig. 3. In order to get a measure of the correspondence between 
the experimental results and the calculated light-cone wave functions, we fit the results with 
a linear combination of squares of the two wave functions. This assumes an incoherent 
combination of the two wave functions and that the evolution of the CZ function is slow (as 
stated in [7]). It is hard to justify these two assumptions because it is hard to make model 
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independent evolutions and to know the phase between the two amplitudes. We therefore 
regard this fit as a qualitative indication of how well each function describes the data. We 
use results of the simulated wave functions (squared) after they were subjected to effects of 
experimental acceptance (Fig. 2, right). 




FIG. 3. The x distribution of diffractive di-jets from the platinum target for 1.25 < kt < 1.5 
GeV/c (left) and for 1.5 < k t < 2.5 GeV/c (right). The solid line is a fit to a combination of the 
asymptotic and CZ wave functions. The dashed line shows the contribution from the asymptotic 
function and the dotted line that of the CZ function. 
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TABLES 



k t bin 
GeV/c 




^stat 

G>as 


^sys 

"as 




O'cz 


Astat 

a cz 


Ucz 




1.25 - 1.5 


0.64 


±0.12 


±0.07 -0.01 


±0.14 -0.12 


0.36 


T0.12 


-0.07 ±0.01 


-0.14 ±0.12 


1.5 - 2.5 


1.00 


±0.10 


±0.00 -0.10 


±0.10 -0.14 


0.00 


TO. 10 


-0.00 ±0.10 


-0.10 ±0.14 



TABLE I. Contributions from the asymptotic (a as ) and CZ (a cz ) wave functions to a fit to the 
data. 







The measured x distributions are shown in Fig. 3 with the combinations resulting from 
the above fits superimposed on the data. The individual contributions from each wave 
function are shown as well. In addition to the statistical errors of the fit, we considered 
systematic uncertainties originating from the background subtraction, from the quality of 
the jets and their identification and from using discrete-mass simulations. The dominant 
contribution comes from the quality of the jets in the low k t region and from using discrete- 
mass MC in the high k t region. The results of the fits are given in Table I in terms of the 
coefficients a as and a cz representing the contributions of the asymptotic and CZ functions, 
respectively. The total errors are obtained by adding the statistical and systematic errors 
in quadrature. 

The values of x 2 /dof were 1.5 and 1.0 for the low and high k t bins, respectively. The re- 
sults for the higher k t window show clearly that the asymptotic wave function describes the 
data very well. Because of the dominance of the asymptotic wave function, this conclusion 
does not depend on the assumptions made in fitting the data to a combination of the two 
functions. The distribution in the lower window is consistent with a significant contribution 
from the Chernyak-Zhitnitsky wave function. However, as stated above, if neither function is 
dominant this can only indicate that at low kt neither function alone describes the data well. 
The requirement that k t > 1.5 GeV/c can be translated (Eqn. 5) to Q 2 > ~ 10 (GeV/c) 2 . 
This shows that for these Q 2 values, the perturbative QCD approach that led to construction 
of the asymptotic wave function is reasonable. 

The kt dependence of diffractive di-jets depends on the quark distribution amplitude. It 
was calculated recently by Frankfurt, Miller, and Strikman [14]. They show that the most 
important terms are those in which the \qq) component of the pion interacts with two gluons 
emitted by the target. The predicted k t dependence can be seen from the cross section for 
this process [14]: 



da 

dk? 



oc 



a s (k 2 )x Bj G(x Bj , k 2 t ) 



d 2 



dk 2 



ip(x, k 2 ) 



(6) 



where ip is the light-cone wave function. Evaluation of the light-cone wave function at 
large kt as due to one gluon exchange gives ip ~ p with a slow function of kt (e.g. the 
asymptotic function). Given the weak a s {k 2 ) dependence, differentiating and squaring gives 



k t 8 . Since a s (k 2 )G(x B j, k 2 ) ~ k]^ 2 [16] the expected dependence of the cross section on 



k t is: 



da 
dkt 



k-\ 



(7) 



This prediction can be compared with the data. We use the MC simulations discussed 
above for di-jets having masses of 4, 5, and 6 GeV/c 2 and the asymptotic wave function 
to correct for the experiment acceptance of the k t distribution. The corrected results are 
shown in Fig. 4(a). Superimposed on the data is a power- law fit k,™ for k t > 1.25 GeV/c. 
We find n = — 9.2 ±0A(stat) ±0.3(sys) with x 2 /dof = 1.0. This slope is significantly larger 
than expected. We note, however, that above k t ~ 1.8 GeV/c the slope changes (although 



in 



the statistical precision there is poor). A power law fit to this region (Fig. 4(b)) results in 
n = —6.5 ± 2.0 with x 2 /dof = 0.8, consistent with the predictions. This would support the 
evaluation of the light-cone wave function at large k t as due to one gluon exchange. 



The steep /^-dependence in lower k t region may be interpreted [14] as a manifestation of 
non-perturbative effects. We try the non-perturbative Gaussian function: ip ~ e _/3fct [17]. 
When inserted in Eqn. 6 and using a s (k 2 )G(x B j-, k 2 ) ~ k^ 2 , we get: 

^ = C(k 2 - 2pkt + (5 2 k\)e~ 2 ^ (8) 

where C is a normalization factor. In Fig. 4(b) we show a fit to this function in the low 
k t range yielding = 1.78 ± 0.05(stat) ± O.l(sys) with x 2 /dof = 1.1. To the best of 
our knowledge, there is no previous report of a direct measurement of (3. Model-dependent 
values in the range of 0.9 - 4.0 were used in calculations of the 7r — 7 transition form factors 
[17]. The present results may indicate that non-perturbative effects are noticeable up to 
k t ~ 1.5 GeV/c, as is the case for the light-cone wave function which becomes dominated 
by the asymptotic function only for larger k t values. 

In summary, we have presented results of direct measurements of the valence quark and 
antiquark momentum distributions in the pion. They show that above k t ~ 1.5 GeV/c 
(Q 2 ~ 10 GeV/c 2 ) the asymptotic distribution amplitude calculated using perturbative 
QCD is applicable. The measured kt distribution in this region is also consistent with this 
conclusion. In the lower k t region there may be other contributions, such as from the CZ 
wave function or other nonperturbative effects. 
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FIG. 4. Comparison of the k t distribution of acceptance-corrected data with fits to cross section 
dependence (a) according to a power law, (b) based on a nonperturbative Gaussian wave function 
for low kt and a power law, as expected from perturbative calculations, for high k t . 
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